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It is shown that there is a double counting in the standard model 
of nn mixing in the medium, resulting in full cancellation of leading 
terms. The direct calculation of nn transition followed by annihilation 
. is performed. The lower limit for the free-space nn oscillation time 

is T n n ~ T n n > 10 16 y, where T nfi is the lifetime of neutron bound 
in a nucleus. This limit exceeds the previous one by 16 orders of 
j_j ■ magnitude. 
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Any information on the occurrence of nn oscillation [1,2] is important in 
order to discriminate among various grand unified theories. The most di- 
rect limit on the free-space nn oscillation time r nn is obtained using free 
neutrons: r nn > 10 7 s[3]. Alternatively, a limit can be extracted from the 
nuclear annihilation lifetime measured in proton-decay type experiments: 
T nfl > 10 s s ~ ly (see, for example, Ref. [4,5]). The process under consid- 
eration is (nucleus) — > (n — nucleus) — > (annihilation products). The 
calculations were based on the potential model of nn mixing in the medium 
[4] or on nonrelativistic diagram technique [5]. They predict the drastical 
suppression of nn transition in nuclei. However, nn conversion comes from 
the exchange of Higgs with rrtu > 10 5 GeV, so from the point of view of the 
microscopic theory (dynamic nn conversion [2], annihilation) the reasons for 
the suppression are not clear. 

In this paper it is shown that the models used previously [4,5] are too 
crude or, more categorically, inapplicable to the problem under study. In 
particular, the potential model does not correspond to annihilation prod- 
ucts in the final state. We perform the direct calculation of the process 
(nucleus) — > (n — nucleus) — > (annihilation products). The interaction 
Hamiltonian is taken in the general form. Then the basic part of calculation 
is a model-independent. However, the amplitude in this case is singular. For 
solving the problem the approach with finite time interval [6] is used. 

l.In the standard approach (later on reffered to as potential model) the 
nn transitions in the medium are described by Schrodinger equations 

(id t + V 2 /2m — U n )n(x) = eii(x), 

(id t + V 2 /2m-U n )n(x) = en(x). (1) 

Here e = l/r n ^ is a small parameter [4]; U n and Un are the self-consistent 
neutron potential and the n-nucleus optical potential respectively. For U n = 
const, and U n = const, in the lowest order in e the probability of the process 
is 

W pot (t) = 1- | U u (t) \ 2 = 2ImT u (t), 
T u (t) = i(e/5U) 2 [l - iSUt - exp(-i5Ut)}, (2) 

where U and T are the evolution operator and T-operator respectively; U = 
1 + iT and 

6U = Un-U n = ReUn - iT/2 - U n . (3) 
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Here T ~ lOOMeV is the annihilation width of n-nucleus state. 

What is described by W pot (t)7 Let us take the imaginary part of Eq.(16) 
ofRef.[6] 

2ImT u {t) = e 2 t 2 - e 2 [ dt a f dt^ImT^r), (4) 

Jo Jo 

where r = t a — tp. Here T(t) and T n (r) are the T-operators of the whole 
process and n-nucleus interaction respectively. For T n -operator \ i>=\0n p > 
and < / | are the n-nucleus and annihilation products respectively. For the 
whole T-operator | % >—\ 0n p > is the nucleus, the physical meaning of the final 
states will be cleared up later. The probability conservation J2f I Ufi \ 2 = 1 
gives us 

= ]T | T fi | 2 + | Ta | 2 , 
E I T fi (t) \ 2 = Wit). (5) 

The probability of the process W{t) will be defined below; | | 2 ~ e 2 , 
whereas | T u | 2 ~ e 4 (see Eqs.(2) and (15)). So for the l.h.s. of Eq.(4) 
2ImT u (t) = W(t), that was taken into account in (2). For the T-matrix of 
T n (r )-operator Eq.(5) has the form 

2Mf (r) = E I T%{r) | 2 + | 7*(r) | 2 , (6) 

T%=< 0n p | T n |0n p >, =< / | T™ | 0n p >. The n-nucleus interaction is 
nonperturbative process and | T£ | 2 ~ I ^/i | 2 - Now Eq.(4) is 



= t 2 t 2 -e 2 ^ dt a ^ dtp |<0n p | T"(r) |0n p >| 2 -e 2 ^ * dt a ^ "dt^E l</l T "( r ) I°™p>| 2 • 

^ 

(7) 

Let us calculate and in the framework of potential model. The 
wave function of initial state is described by the equation 

i w = H0, (8, 

Ho = — V 2 /2m + C/ n . At the moment t — the interaction 5C/ is turned on. 
We have 

i— = (^o + ^)*, (9) 
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^(0) = $(0). The projection to the initial state and T- matrix at t — r are 

< $ | q? >= u*(t) = exp(-i8UT), (10) 
7^(r)=i[l-expH^r)], 
E I T fi( r ) \ 2 = I" I | 2 = 1 - e"^ = W n (r), (11) 

where W Ti (t) is the n-nucleus decay probability. Note that Y corresponds to 
all n-nucleus interactions followed by annihilation. However, the main con- 
tribution comes from the annihilation without rescattering of n[5], because 
<7 ann > 2a sc . Substituting these expressions in (7), one obtains the potential 
model result (2). 

Therefore, the approach with the finite time interval was verified by the 
example of exactly solvable potential model. It is involved in Eq.(4) as a 
special case. Solving Eqs.(l) by method of Green functions we will obtain 
the same results. We have started from Eq.(4) only for verification of the 
finite time approach. 

Let us return to Eq.(7). It is at least unclear. (l)The first term is 
free-space nn transition probability. The matrix elements T£ and Tj^ de- 
scribe transitions (n — nucleus) — > (n — nucleus) and (n — nucleus) — > 
(annihilation products) respectively. So the first and the second terms cor- 
responds to n-nucleus in the final states. However, in the experiment only 
the annihilation products are detected (n-nucleus is unobservable) and the 
result should be expressed only in the terms of T^,^ Moreover, n-nucleus 
decays into final states identical with the states given by the third term. This 
suggests that the potential model contains the double counting. Expression 
1— | Uu | 2 from Eq.(2) describes the inclusive decay of initial state and so the 
nn transition with n-nucleus in the final state is also included in W pot , unless 
additional limits are imposed. To exclude the double counting the annihila- 
tion products in the final state should be fixed. (2) Let us | SUt |<C 1. (This 
is the case in some other problems.) When T = 0, the third term equals 
to zero. When r ^ 0, the contribution of the third term is negative and 
dW/dT < 0, whereas the opening of the new channel (annihilation) should 
increase W. 

How big is the probable error? The contributions of the second and third 
terms are: x 2 = — e 2 t 2 /2 + F 2 , x 3 = —e 2 t 2 /2 + F 3 . The functions F 2i3 contain 
the terms proportional to t and exp(— iSUt). So the e 2 t 2 term produced by 
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the third term is fully canceled. This is a consequence of double counting. 
Therein lies the reason of the discrepancy between our result and the result 
of the potential model. 

As noted in [6], Eqs.(ll) and (2) can also be obtained by means of micro- 
scopic variant of the potential model (zero angle rescattering diagrams of n). 
In this case the Hamiltonian of n-medium interaction is H = SU. The same 
calculation was repeated by Dover et al. [4]. They substitute H = —iT/2 in 
(4) and obtaine (2). On the basis of this and only this they refute the result 
of Ref. [6]. In other words they refute our limit because it differs from the 
prediction of the potential model (H = —iT/2[A]). 

What is wanted is J2f=£i I | 2 , where </| is the annihilation products. 
It is connected with the diagonal matrix element by Eq.(5): 



Calculation of Ta is determined by r.h.s. of Eq.(12): the cut corresponding 
to Tu must contain only annihilation products, that is not in accordance 
with Eq.(7). It includes redundent states / = (n — nucleus) forbidden by the 
unitarity condition. The relation (12) is not fulfilled. Also the eigenfunctions 
of H + SU do not form the complete orthogonal set. Due to this fact the n- 
nucleus (described by U n ) also can not appear in Eq.(12) as the intermediate 
state. So the model (1) is inapplicable in our case because it leads automat- 
ically to incorrect matrix element T ti . Elimination of redundent trajectories 
from Tjj means the direct calculation of T^. 

2. In Ref. [6] the first and the third terms were taken into account. The 
second one was omitted. The first term corresponds to low density limit 
and is meaningfull for nn transitions in the gas. This scheme is not quite 
correct here. In this paper we perform the direct calculation of the process 
{nucleus) — > (n — nucleus) — > {annihilation products). We have 



2ImT ii = 'ET%T fi . 



(12) 



t 



f^ 1 dt k /** dtpH(t 1 )...H(t k )H nfi (t /3 ) 



< f | U{t,0)-I | 0n p >= iT fi (t) = EH)" +1 <f\ dti- 
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fc=i 



(13) 



where 



H{i) = {all n — medium interactions) — U, 




(14) 
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H + H nn = Hi. Here | On p > is the state of the medium containing the neu- 
tron with 4- momenta p = (p^/2m + U n , p n ), < / | represents the annihilation 
products; H nn is the oscillation Hamiltonian[4]. In the case of the formulation 
of the S'-matrix problem (t, 0) — > (oo, — oo) Eq.(13) in the momentum repre- 
sentation includes the singular propagator G = l/(e n — p^ t /2m — U n ) ~ 1/0. 
Taking into account that H nfl \ 0n p >— e \ 0n p >, we change the order of 
integration and obtain 

T fi (t) = -e f 'dtpiT*(t-t p ), 

iT%{r) = J2(-*) k / dh... / dt k <f\ H(t 1 )...H(t k ) |0n p >, (15) 
fc=i Jt Jt P 

where | 0n p > is the state of the medium containing the n with 4-momenta p; 
t = t — tp. The 4-momenta of n and n are equal. is an exact amplitude 
of n-nucleus decay. It includes all the n-nucleus interactions followed by 
annihilation. Expression for Tfi(t) was obtained in perfect analogy to (4) 
that can be considered test for Eq.(15). 

The 2-step process was reduced to the annihilation decay of n-nucleus. 
(The slightly different method is the separation of the antineutron Green 
function [6].) It is seen from (13) and (15) that both pre- and post- nn 
conversion spatial wave functions of the system coincide: 

\0n p > sp =\0n p > sp . (16) 

n appears in the state with SU = 0. We would like to stress that in the 
potential model (1) the picture of n-nucleus formation is exactely the same: 
in Eq.(4) for T u =<0n p | T \ 0n p > and T% =<0n p | T n \ 0n p > condition (16) 
was fulfilled. Hereafter, the potential model of the n- medium interaction 
(block T n ) was used and W pot was reproduced, which confirms the picture 
of n-nucleus formation given above. Solving Eqs.(l) by method of Green 
functions we will obtain the same results, including (16). The equality of 
vectors of state (16) is also evident from the continuity of solution of Eqs.(l). 

In both models the first stage of the process (nn conversion) is described 
identically. The basic difference centers on the next stage - annihilation. In 
the potential model T£ is calculated (as a result the self-energy part E = 5U 
appears) and is used in Eq.(7), which is wrong. We calculate start- 
ing from the same point (16). The result will be expressed through T (see 



6 



Eqs.(19),(ll)), but not through SU = ReSU — iT/2, as is usually the case in 
decay calculations. The standard ^[/-dependence is manifested in scattering 
problems, when the diagonal matrix element Ta in r.h.s. of Eq.(12) should 
be taken into account. It corresponds to the observarble process - zero angle 
scattering of the incident particle. 

The characteristic annihilation time of n is: A = 1/T ~ 10~ 23 s. When 
t ^> A, T^(r) reaches its asimptotic value T^: 

T%{r » A) = T"(oc) = T% = const. (17) 

The expressions of this type are the basis for all S'-matrix calculations. (Mea- 
surement of any process corresponds to some interval r. So it is necessary to 
calculate U(t). The replacement U(r) — > S(oo) is equivalent to (17).) Let 
us take t ^> A. From (15) and (17) we have 

T fi (t) = -ie[ dt p T%{t - t p ) + f dt p T%(t - t p )\ ~ -ietT%. (18) 

The contribution of the second term is negligible since | T^(r) | 2 < 1. The 
probability of the whole process is 

W{t) = E | T fi (t) | 2 ~ e 2 t 2 I T?M \ 2 - e ^ 2 ' (19) 

where Eq.(ll) has been taken into account. The value e 2 t 2 = t 2 /r 2 n is the 
free-space nn transition probability. Due to the annihilation channel nn 
conversion is practically unaffected by the medium. So r nn ~ T nn , where T nfl 
is the oscillation time of neutron bound in a nucleus. In order to find the 
limit for r nn from experimental data on nuclear stability, the distribution (19) 
should be used (but not the exponential decay law!). Let us N n ,T ,ei and 
6 are the total number of neutrons under observation, the observation time, 
the overal n n detection efficiency and the average number of observable 
n — > n events respecrively. From the inequality 

iVn(T /r nfi ) 2 ( ei /#) < 1 (20) 

one obtains r nn > 10 w y, where the values T = 1.3y, N n = 2.4- 10 32 , e x = 0.33 
and 6 = 2.3 [7] were used. 

Our previous result [6] is different from (19) only by a factor of 2. How- 
ever, in Ref. [6] we used the limit T nn > 4.3 • 10 31 y [7] deduced from the 
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experimental data by means of exponential decay law which does not agree 
with (19). 

3. Let us return to the reason of enormous quantitative disagreement be- 
tween our result and the potential model one. The strong sensitivity of the 
results should be expected. Really, in the momentum representation the S- 
matrix amplitude M s , corresponding to nn transition followed by annihilation 
(see Eq.(14)) diverges 

M s = e -J M ~ -, (21) 

e n -p 2 n /2m-U n 0' V ; 

where M is the annihilation amplitude. These are infrared singularities con- 
ditioned by zero momentum transfer in the e-vertex. It is can be seen that 
M s ~ 1/0 for any bound state wave function of a neutron (i.e., for any nu- 
clear model). On the other hand in the potential model the energy is not 
conserved and becomes complex: Ma — > Ma + SU (Ma is the nuclear mass). 
In this case we have instead of (21) 

M pot = e -L — M = e^-M. (22) 

This is a potential model amplitude. Really, the process width is T pot — J d$ \ 
M pot | 2 /2M4 = e 2 r/ I SU | 2 = W pot /t, that coincides with (2) when | SUt |> 
1. It is seen that: (1) There is a double counting in M and G with respect to 
H . M po t does not agree with Eq.(14) as well. (2) SU = is the singular point 
and due to zero momentum transfer q = in the vertex corresponding to 
H nn we are in this point. So the result is extremely sensitive to SU. (Usually, 
in the reactions and decays the momentum transferred is q 7^ 0. In this case 
^[/-dependence of G is masked by q: G^ 1 = (e n —qo) — (p n —q) 2 /2m—U n —SU. 
We deal with the 2-tail and q — 0.) 

Comparing (21) with (18) one sees that in principle the limit SU — > 
corresponds to the replacement 

1/SU -> t. (23) 

Certainly, we do not set SU = 0. U n is not introduced at all. In the calcu- 
lation of Eq.(13) the multiplier t (see Eq.(18)) arises automatically instead 
of 1/SU in the potential model, or 1/Aq in the case g ^ 0. When g ^ 
in the e-vertex, Eq.(13) leads to usual ^-matrix result (see below). Formal 
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reason for the differences in the results is the full cancellation of the terms 
~ t 2 in Eq.(7). Erroneous structure of (7) is caused by nonperturbative and 
2-step nature of the process. The fact that q = extremely increases the 
disagreement. 

4. An additional comment is necessary regarding t-dependence of the whole 
process probability W(t). Eq.(19) has been obtained in the lowest order in 
e. The exact distribution W pr (t) which accounts for the all orders in e is 
unknown. However, W is the first term of the expansion of W pr and we can 
restrict ourselves to the lowest order W pr = W, as it is usully the case for 
rare decays. W po t is also calculated in the lowest order in e. 

The protons must be in very early stage of the decay process. Thus the 
realistic possibility is considered [8-11] that the proton has not yet entered 
the exponential stage of its decay but is, instead, subject to non-exponential 
behavior which is rigorously demanded by quantum theory for sufficiently 
early times. At first sight, since T nfl > 10 16 y for nn-mixing in a nuclear, 
the non-exponential behavior should be expected too. In fact, there is one 
more problem: we deal with the two-step process. When trying to calculate 
M s and T s in the framework of standard S-matrix theory we get T s ~ 1/0. 
So the decay law exp(— T s t) is unrelevant and it is necessary to deduce the 
distribution W(t) as it was done above. 

We have to mention the main points of Krivoruchenko's preprint [12]. 
(l)The nn transition followed by annihilation (two-step nuclear decay) and 
motion of particle in the classical field are two different problems. Describ- 
ing the first one by Eqs.(l) we understand that this is an effective procedure. 
From formal standpoint in the first and the second cases the potentials are 
complex and real respectively. Unfortunately, sometimes the literal analogy 
between these problems is drawn[12]. (2)The initial Eq.(ll) of Ref.[12] must 
describe the nn transition followed by annihilation. However, the l.h.s. of 
Eq.(ll) is free of n-nucleus interaction at all. The r.h.s. contains annihila- 
tion width T (we stress this point) and coincides with the potential model 
result. We also would be glad to get the result without calculations, but 
some difficulties emerge in reaching this goal. 

The interaction responsible for the nn conversion is ultra-weak. There- 
fore, the n-nucleus interaction in the initial state should be taken into account 
exactly. The neutron line entering into the nn transition vertex should be 
the wave function of the bound state (see Eq.(8)), but not the propagator, 
as in the model based on diagram technique [5, 13]. As a result, in this model 
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the nn transition is possible only between the acts of interactions of oscillat- 
ing particle and a nucleus. These interactions lead to total suppresion of nn 
conversion, that is incorrect. This can be understood using the analogy with 
(3 decay and taking into account that m H > 10 5 GeV. 

Some additional remarks to paper[13] are necessary. (l)The picture de- 
scribed in Sec.l is valid only for simple interaction operator. We deal with 
products of operators (see Eq.(13) of this paper). (2) In Sec. 2 it is claimed 
that the amplitude should be singular (!) at B — > 0, where B is the binding 
energy. Accordingly, as B — > 0, the amplitude obtained | A\ | 2 ~ 1/0. In fact, 
they are usual infrared singularities mentioned above, wich must be avoided. 
The correct model should reproduce the law dencity limit W(t) = e 2 t 2 . 
(3)The cut corresponding to the diagonal matrix element (18) is completely 
free of annihilation products. 

We try to calculate the process amplitude starting from (14). The S- 
matrix theory gives (21). The approach with finite time interval is infrared- 
free. Its verification for the diagrams with q = was made above by the 
example of potential model. For nonsingular diagrams the test is obvious. Let 
us q 7^ in the e- vertex. The appropriate calculations with finite time interval 
(adiabatic hypothesis should be used) give the ^-matrix result (we stress this 
fact because it means the verification of the approach): T/« = ie'(l/ Aq)Tjl, 
where TJ\ is the S'-matrix amplitude of annihilation of virtual n with 4- 
momenta k = p — q. Comparing with (18) one sees that limit Aq — > 
corresponds to the replacement l/Aq — > t,Tj\ — > (compare with (23)). 
The similar problem for matrix element Ta was solved in Ref.[14]. 

The main results of this paper are given in the abstract. In the next 
paper the following statements will be proved. (1) All the results are true 
for any nuclear model. (2) The contribution of the corrections is negligible. 
(3) Further investigation and verification of the approach will be presented 
as well. In our opinion, it makes sense to look at some other problems on 
oscillation of particles in a medium from the standpoint given above. 
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